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ABSTRACT 

The double pulsar PSR J0737-3039A/B displays short, 30 s eclipses that arise around conjunction when the 
radio waves emitted by pulsar A are absorbed as they propagate through the magnetosphere of its companion 
pulsar B. These eclipses offer a unique opportunity to probe directly the magnetospheric structure and the 
plasma properties of pulsar B. We have performed a comprehensive analysis of the eclipse phenomenology 
using multi-frequency radio observations obtained with the Green Bank Telescope. W e have characterized 
the periodic flux modulations previously discovered at 820 MHz by iMcLaughlin et all and investigated the 
radio frequency dependence of the duration and depth of the eclipses. Based on their weak radio frequency 
evolution, we conclude that the plasma in pulsar B's magnetosphere requires a large multiplicity factor (~ 10 5 ). 
We also found that, as expected, flux modulations are present at all radio frequencies in which eclipses can be 
detected. Their complex behavior is consi stent with the confinem ent of the absorbing plasma in the dipolar 
magnetic field of pulsar B as suggested by iLyutikov & Thompsonl and such a geometric connection explains 
that the observed periodicity is harmonically related to pulsar B's spin frequency. We observe that the eclipses 
require a sharp transition region beyond which the plasma density drops off abruptly. Such a region defines a 
plasmasphere which would be well inside the magnetospheric boundary of an undisturbed pulsar. It is also two 
times smaller than the expected standoff radius calculated using the balance of the wind pressure from pulsar 
A and the nominally estimated magnetic pressure of pulsar B. 

Subject headings: binaries: eclipsing — pulsars: individual (PSR J0737-3039A/B) 



1. INTRODUCTION 

The double pulsar PSR J0737-3039A/B is one of the most 
celebrated systems in astrophysics. On the one hand, the fact 
that it comprises two neutron stars, both visible as radio pul- 
sars, combined with being the most relativistic pulsar binary, 
allows for the measurement of post -Keplerian paramete rs at 
an unprecedented level of accuracy dKramer et al.ll2006h . Its 
short coalescence timescale significantly increases the esti- 
mates of the expe cted double neutron star merging rate in 
galaxies like ours dBurgav et al.l 120031) . the timing provides 
the most prec ise test of general r elativity in the strong-field 
regime so far (iKramer et al.ll2006h and, more recently, a high- 
precision measurement of the relativistic spin precession of 
pulsar B yielded the first theory-independent constraint on 
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spin-o rbit coupling in the strong-field regime dBreton et alj 
l2008h . 

On the other hand, the double pulsar offers a unique labora- 
tory to study a variety of pulsar phenomena. This is related to 
the fact that pulsar "A", the 23 ms spin period recycled pulsar, 
has an energetic wind roughly 3625 times stronger than that 
of its 2.8 s companion, "B". As a consequence, the dynamic 
pressure from pulsar A's wind confines and distorts the mag- 
netosphere of pulsar B, whose radio pulsations are difficult 
to detec t throughout mos t of the orbit except in two bright 
regions dLvne et all 120041) . It was proposed that pulsar A's 
wind changes the direction of B' s radio emission in different 
parts of the orbit (Lyutikov 2005), but the exact nature of the 
phenomeno n still remains uncertain despite pheno menologi- 
cal studies dBurgav et alJl2005t iPerera et alJl2010h . Pulsar B 
also shows sub-pulse drifting at the spin fre quency of pulsar A 
in one of the bright orbital phase windows dMcLaughlin et all 
l2004bl) . demonstrating a direct modulation of B's emission 
by the non-axisy mmetric dynamic p ressure of A's Poynting- 
dominated wind dSpitkovskyl 120061) . It might even be pos- 
sible to use this e ffect as an extra tool to time the system 
dFreire et alj|2009l) . 

Arguably, the most impressive display of the interaction be- 
tween the two pulsars are the radio eclipses of pulsar A when 
it passes behind pulsar B, thanks to the fortunate nearly edge- 
on viewing geometry. As opposed to eclipses in other types of 
pulsar binaries, such as those resu lting from the stellar w ind of 
a main sequence companion (e.g. Uohnston et alii 19941) or the 
abla ted material from the surface of a low-mas s companion 
(e.g. lThorsett & Nicell99UIStappers et"aT]|2001l) . those of the 
double pulsar have a magnetospheric origin. The duration of 
the eclipses implies a sliced transverse size of 18, 100 km if the 
absorbing region lies around pulsar B at a distance of 2.9 lt- 
s from pulsar A dKaspi et alj |2004): this is well within pul- 
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sar B's inferred light cylinder radius, which is 135,000km, as 
well as the estimated 40,000 km magnetosphere radius in the 
presence of pulsar A's wind. This first occurrence of eclipses 
in a double neutron star system provides an exceptional way 
of probing the magnetosphere of a pulsar. 

One of the most outstanding problems in understanding 
pulsars is the generation of coherent radio emission across 
a wide range of rotational periods and magnetic fields. The 
physics behind the emission is intimately linked to the prop- 
erties of their magnetospheres, in particular their structure and 
the particle density, which have yet to be well constrained 
observationally and theoreti cally. A realistic magneto sphere 
may deviate from a classical lGoldreich & Julian! (11969b dipo- 
lar magnetosphere in having a particle density exceeding the 
so-calle d Goldreich-Julian density, exc ept very close to its 
surface dRuderman & Sutherland! fl 975). Theoretical work 
shows that the creation of additional pair plasma in the mag- 
netosphere is a key ingredient to power the radio emission but 
progr ess still has to be made to reconcile theory and observa- 
tions dAronsll2000l) . In the case of the double pulsar, we have 
for the first time a powerful source of radio light positioned 
behind, albeit externa lly modified, a pulsar magnet osphere. 

As pointed out by iRafikov & Goldreichl ( 120051) . resonant 
cyclotron absorption should be the dominant source of extinc- 
tion of radio waves in an isolated pulsar magnetosphere but 
this mechanism would only yield optical depths in the range 
of unity if the multiplicity factor is the order of 100. One 
way to reach such a value is if the particle number density 
and trapping efficiency can be increased via the interaction 
with pulsar A. As a consequence, though, particles will be 
heated to relativistic velocities and resonant synchrotron ab- 
sorption will inevitably become the main source of opacity. 
While resonant synchrotron absorption appears to be a viable 
eclipse mechanism, a number of questions still remain. How 
large can the electron multiplicity factor become? What is the 
radial density profile of the electrons? Through which mech- 
anism is the heating resulting from the interaction with pulsar 
A operating? 

In their model, Rafikov & Goldreich ( 2005, hereafter 
RG05) propose that if the illumination from pulsar A's ra- 
dio beam can pump sufficient perpendicular momentum in the 
particles of pulsar B's magnetosphere via cyclotron absorp- 
tion that the trapping efficiency can increase significantly. In 
turn, the particle density would be able to grow and reach a 
steady state. Their model predicts a multiplicity factor ^100 
at a radius 10,000 km for fiducial neutron star parameters. It 
also predicts an eclipse duration scaling as v~ x l l , where v is 
the observed radio frequency. 

In c ontrast, the eclipse model from iLvutikov & Thompson] 
d2005l hereafter LT05I), suggests that the multiplicity factor 
is much larger, ~ 10 5 . They believe it can be provided by 
torsional Alfven waves resulting from the interaction of pul- 
sar B's magnetosphere with pulsar A's relativistic wind. A 
particular feature of their model is that particles should be 
trapped in the closed field lines of the magnetosphere and 
have a roughly constant density. To simplify the model fur- 
ther, they suggest that the density distribution can be assumed 
constant up to the field lines closing at some characteristic 
radius, beyond which trapping becomes inefficient and the 
density drops to negligible values. Alternatively, they also 
suggest a different version of their model, hereafter LT05II, 
that retains the overall geometrical properties but in which 
the absorbing electrons are supported by an ion cloud. In 
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this model, the electron number density is constant along each 
field line and proportional to the electron temperature, which 
is itself a function of the radius (suggested to be a power law, 
T e {r) <x r 1 ). 

In this paper, we revisit some of the phenomenological 
properties of the double pulsar eclipses in light of the progress 
made in our understanding of the phenomenon using the 
LT05I model. Indeed, this model, which relies on a mag- 
netic dipole misaligned with respect to the spin axis, has 
be en used successfully to reproduce the eclipse light curves 
by iBreton et al.l (120081) and to explain the B's pulse profile 
changes bv lPerera et alj (120101) . After introducing our obser- 
vations and data reduction in Section [2] we present in Sec- 
tion[3]a multi-radio-frequency analysis of the peculiar eclipse 
properties including flux modulations, their periodicity, du- 
ration and depth. Whereas the periodic behavior of the flux 
modulations is tightly connected to the geometry of pulsar B 
with respect to the orbit, we show that the eclipse duration 
and depth can be used to probe the properties of the plasma 
confined in the magnetosphere. In Section [5] we summarize 
our work and provide some concluding remarks. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the double pulsar regularly since 2003 De- 
cember as part of a multi-purpose monitoring campaign aimed 
primarily at the timing of the pulsars, but also at the investiga- 
tion of the different phenomena displayed by this system such 
as the eclipses of pulsar A, the orbital modulation of pulsar 
B's flux, and its subpulse drifting. Data used for the analysis 
presented in this paper were acquired at the Green Bank Tele- 
scope (GBT), in West Virginia, with the SPIGOT and BCPM 
backend instruments dKaplan et al.H2005b iBackeret alii 19971) 
between 2003 December and 2007 November (proposal codes 
GBT03C-041, GBT04B-026, GBT05B-042, GBT06B-018 
and GBT07B-029). For the multi-frequency analysis re- 
ported in the next section, we use only the observations ob- 
tained during the first year (MJDs 52984-53379) because 
they cover a large range of radio frequencies while also span- 
ning a short enough time that the eclipse profile did not 
change significantly due to the effect of relativistic precession 
of pulsar B's spin axis. Relevant details about the observa- 
tions and setup are presented in Tables Q] and [2] 

We performed the initial data reduction using the pul- 
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sar analysispackages PRESTcQ dRansom et alj 120021) and 
SIGPROC0 First, we dedispersed the data (DM = 48.920 
cm" 3 pc) to correct for the frequency-dependent travel time in 
the ionized interstellar medium. Then we generated folded 
data products using the predicted spin periods of pulsars A 
and B. F or these two s t eps, w e used the timing solution pre- 
sented in lKramer et al.1 (120061) . 

We generated a pulsed flux time series for each segment 
of observations cont aining an ecl i pse of pulsar A in the same 
way as presented in iBreton et al.l 120081) . We generated high 
signal-to-noise (S/N) pulse templates of pulsar A by integrat- 
ing its flux over each observation. Then we calculated the 
relative pulsed flux density of pulsar A by fitting the corre- 
sponding pulse profile templates to folded intervals consist- 
ing of the sum of four individual pulses of pulsar A (~91 ms). 
This choice provides a good balance between S/N and time 
resolution of the eclipse features. Note that as opposed 
to SIGPROC, which folds at an integer number of pulses, 
PRESTO's folding algorithm is based on the data size and 
hence despite the fact that the exact number of pulses is gener- 
ally four, there are a handful of data points consisting of three 
pulses only. We renormalized these data points to correct for 
their shorter integration time. Finally, we normalized the time 
series such that the flux average in the out-of-eclipse region is 
unity. 

Radio frequency interference is an issue that limits the qual- 
ity of our data. We observe that, with similar setup and ob- 
serving conditions, the S/N can vary substantially from one 
observation to another depending on the level of interference. 
Some of it can be mitigated if it is well localized in the radio 
spectrum but our determination of the pulsed flux amplitude is 
mainly affected by short interference "bursts", which are im- 
possible to eliminate, especially because we deal with short 
folding intervals. Generally speaking, GBT observations at 
820 MHz are relatively free of interference and this motivates 
why data are often recorded at this frequency . 

Using the timing solution of iKramer et alJ (120061) . we cal- 
culated the orbital phase at each data point of our flux time 
series. Throughout this paper, we refer to the orbital phases 
using the superior conjunction of pulsar A (i.e. the eclipse 
mid-point) as zero-point. We also empirically determined the 
spin phases of pulsar B associated with the data points of pul- 
sar A's time series using the method described in Section 1 of 




'■ PRESTO is freely available at http : / /www . cv . nrao . edu/~s ransom/presto 
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Figure 1. Typical individual eclipse light curves of pulsar A observed with 
the GBT at 325, 427, 820, 1400, and 1950 MHz. Each point in the light 
curves was made of four individual pulses from pulsar A, while the vertical 
dashed lines mark each one of pulsar B's pulsations. At the three lowest radio 
frequencies, the modulation behavior is clearly visible whereas the poorer 
S/N at higher frequencies makes it hard to identify unambiguously. Note that 
these observations present higher-than-normal S/N, mainly due to the level 
of radio frequency interference being low. 

the Supporting Online Material of lBreton et al.l (12008b . 

3. ECLIPSE PHENOMENOLOGY 

3.1. Modulations 

What makes the double pulsar eclipses very different from 
any other eclipses ob served to date are t he rapi d flux vari- 
ations discovered by iMcLaug hlin et al. (2004a) in a high 
time resolution analysis of 820 MHz eclipse light curves. 
They found that not only does pulsar A's flux vary on short 
timescales, but that it is also modulated in phase with the ro- 
tation of pulsar B (see Figure [T). This kind of behavior has 
considerable importance for constraining the eclipse origin. 
Early theoretical work to explain the eclipses involved the 
magnetosheath of pulsar B, which would form in the direc- 
tion opposed to pulsar A as a result of th e relativistic shock 
of its wind with pulsar B's magnetosphere dArons et al.H2005t 
lLvutikovll2004l) . This family of models, however, can now be 
excluded as they are unable to reproduce the rapid modulation 
seen in the eclipses. 

3.1.1. Periodicity 
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is freely available at 



Although the connection between the flux variability and 
the rotational phase of pulsar B has been established at 
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Figure 2. Power spectra of the eclipse light curves presented in Figure fj] 
Data were high-pass filtered in order to remove the low-frequency eclipse 
trend. The vertical dashed lines indicate the first four harmonics of pulsar B's 
spin frequency (z^b) at 0.36, 0.72, 1.08, and 1.44 Hz, respectively. 



820 MHz dMcLaughlin et al.ll2004ah . the presence of such be- 
havior at other radio frequencies, although expected, is yet to 
be confirmed. For intrinsic (i.e. spectrum) and technical (i.e. 
instrument bandwidth, sensitivity and radio interference) rea- 
sons, pulsar A at GBT is generally more strongly detected at 
820 MHz than in other frequency bands. Hence, most of the 
data were collected at this frequency, which was also chosen 
for the original analysis bv lMcLaughlin et al.l d2004al) . 

We performed a Fourier analysis on pulsar A's flux 
lightcurves in order to improve the sensitivity to detect mod- 
ulations at other frequencies, as well as to obtain quantitative 
information about this behavior. First, we removed the over- 
all low-frequency trend of the eclipse using a high-pass filter 
consisting of two running means having a window size of 71 
data points (i.e., 0.3° in orbital phase). We found that the 
power spectra of the individual light curves observed at radio 
frequencies 325, 427, 820, 1400, and 1950 MHz show signif- 
icant excess power at the frequency and higher harmonics of 
the rotational period of pulsar B (see Figure |2). These fea- 
tures are easily visible even in data presenting poor S/N and 
for which the modulation is not clearly identifiable by eye. 

The periodic nature of the modulations, harmonically re- 
lated to the spin frequency of pulsar B, is important because 
it suggests a strong connection between the distribution of 
the absorbing material and the geometry of pulsar B. We 
also searched for modulations at radio frequencies higher than 
1950 MHz but the very poor S/N — at 2200 MHz, S/N- 0.25 
per data point in the out-of-eclipse region — prevented us 
from seeing either an eclipse trend or finding any significant 
feature in the power spectrum. Unfortunately, pulsar A be- 
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Figure 3. Dynamic power spectrum analysis of the average light curve con- 
sisting of eight eclipses observed at epoch MJD 54200. In this figure, the 
frequency of the dynamic power spectrum is displayed vertically as a func- 
tion of orbital phase, horizontally. The color intensity indicates the spectral 
energy density, with black being the highest value, and three isocontour levels 
drawn with solid lines. The right panel shows the "standard" power spectrum 
averaged over orbital phases —2° to 2° . Data were high-pass filtered in or- 
der to remove the low-frequency eclipse trend. The horizontal dashed lines 
indicate the first two harmonics of pulsar B's spin frequency (ug) at 0.36 and 
0.72 Hz, respectively. 

comes so faint at higher radio frequencies that in order to ob- 
tain a reasonable S/N, the binning of the sub-fold intervals 
must be longer than the modulation timescale. 

Another insightful way to analyze the flux modulations is 
through the use of a dynamic (also known as windowed) 
Fourier transform. This technique consists of calculating the 
Fourier transform within subsections of a time series and 
repeating the process after translating the window kernel. 
The corresponding power spectra can then be calculated and 
stacked in a two-dimensional plot so that the dynamic infor- 
mation about the time evolution of the power content can 
be visualized. Figure [3] shows such a dynamic power spec- 
trum obtained from the average of eight 820 MHz observa- 
tions observed around MJD 54200. We chose these data be- 
cause they comprise the largest sequence of 820 MHz data 
that we obtained over a short time span. The averaging of 
the ligh t curves was done as described in iMcLaughlin et al] 
d2004al) . That is, because the phase of the modulation varies 
from one eclipse to another, the light curves were shifted by 
up to ±P S pi n .B in order to sum the features coherently. 

Our dynamic power spectrum analysis yields two impor- 
tant results. First, no periodic flux variations are detected out 
of the eclipse region, where the power content is well charac- 
terized by white noise. Despite the fact that our pulsed flux 
time series are made from data folded at pulsar A's spin pe- 
riod, a process which should average out the signature of pul- 
sar B, the fact that modulations occur only during the eclipse 
confirms they are not an artifact of the data reduction. Sec- 
ond, we find that the power content significantly evolves dur- 
ing the eclipse. As seen in Figure [3] we observe that the 
eclipse ingress is initially characterized by modulations hav- 
ing significant power content at the spin frequency (v%) and 
the second harmonic (2^b) of pulsar B. Then, the modulation 
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Figure 4. Folded eclipse profile for the eight eclipses observed at epoch MJD 
54200. Eclipses were combined before being folded at the spin period of 
pulsar B. Lighter regions have higher flux intensity levels. 

switches to a mode dominated by the second harmonic (2j/b) 
around orbital phase -0.5°. It then moves to a mode charac- 
terized by a high harmonic content, which indicates that the 
modulation features are very narrow. Finally, the modulations 
temporarily stop before resuming in the first harmonic (z^b) at 
the egress. 

3.1.2. Folded Light Curve 

To complete the picture, folding the light curve at the spin 
period of pulsar B provides one more way to visualize the 
properties of the eclipses, in particular the dynamic informa- 
tion of the phase of the modulations. Figure|4]shows the same 
eclipse data as in Figure [3]but without the hig h-pass filtering. 
As pointed out by iMcLaughlin et al.1 (l2004al) . we can clearly 
see that the eclipse is longest when the spin phase of pulsar B 
is 0.0 (i.e., when we receive its radio pulse). In a narrow range 
around pulsar B's spin phase 0.25, pulsar A barely disappears. 
We also observe that initially, most of the modulation power 
is due to the absorption feature located around spin phase 0.0, 
which then also appears around spin phase 0.5. Toward the 
center of the eclipse, at orbital phase 0.0, the modulation con- 
sists mainly of narrow emission features at B pulsar phase 
0.25. 

The complex modulation pattern and the asymmetric shape 
of the eclipse arise naturally as a consequence of the mag- 
netic dipole being misaligned with respect to the spin axis in 
the context of the LT05I geometrical model. As seen in Figure 
[5] if we consider the closed field lines reaching maximum ex- 
tent, the magnetosphere of pulsar B resembles a torus which 
presents its largest cross-section at spin phase 0.0, when its ra- 
dio beam is illuminating us. Conversely, the narrowest side of 
the dipole is visible about a quarter of a rotation later/before. 
At the start of the eclipse, the magnetosphere of pulsar B 
intercepts pulsar A's line of sight twice per rotation while 
it remains visible to us when the magnetosphere is viewed 
from the side. Gradually, as the orbital motion brings pul- 
sar A closer along the line of sight of the projected mag- 
netosphere of pulsar B, A becomes hidden most of the time 



and barely emerges in a narrow transparent window when the 
magnetosphere is viewed sideways. The offset of the mag- 
netic axis with respect to the spin axis, ~ 20° dBreton et al.l 
120081) . is such that the magnetosphere of pulsar B will only 
clear out the line of sight at spin phase 0.25, but not at spin 
phase 0.75, which explains the single modulation per rota- 
tion of B in the light curve around conjunction. The large 
ingress/egress asymmetry appears from the fact that the lon- 
gitude of pulsar B's spin axis, about 50° away from our line- 
of-sight vector in 2004, causes the magnetosphere to intercept 
the line of sight only during the egress for the spin phases 0.25 
and 0.75. Based on the geometry inferred from the truncated 
dipole LT05I model, relativistic spin precession of pulsar B 
should bring its spin axis toward our line of sight and the two 
should lie in a plane perpendicular to the orbit around year 
2013. At this time, the eclipse profile should be nearly sym- 
metrical. 

3.2. Duration 

In principle, the eclipse profile not only offers a way of de- 
termining the geometry of pulsar B but also of studying the 
properties of the plasma inside its magnetosphere. For in- 
stance, the radial distribution of the absorbing particles should 
change the amplitude and shape of the eclipse features. How- 
ever, the fact that the optical depth increases very abruptly to 
large values at the eclipse characteristic radius makes it hard 
to probe the inner regions of the magnetosphere since most of 
the differences that arise in the light curve profile happen on 
a timescale and at a flux level that is beyond the observational 
capabilities of the GBT 

On the other hand, the behavior of the eclipse light curves 
as a function of radio frequency mig ht provide an easie r al- 
ternative to learn about the plasma. iKaspi et ail ( 120041) re- 
ported that the eclipse duration is only slightly frequency de- 
pendent in the range 427- 1400 MHz. From the full width 
at half-maximum (FWHM) of the eclipses, which they calcu- 
lated after fitting a Fermi function to each half of the eclipse, 
they found a small, non-zero linear decrease of duration as 
a function of frequency (-4.52(3) x 10~ 7 orbits MHz" 1 , i.e., 
-1.63(1) x 10" 4o MHz _1 ). Since our data set includes ob- 
servations made over a larger range of radio frequencies, 
350-1 950 MHzQ we conducted a similar analysis in order 
to test the predictions of different eclipse models. 

We chose observations made around the same epoch, be- 
tween 2003 December and 2005 January (MJDs 52984- 
53378), in order to reduce effects of the long-term evolution 
of the ecl ipse profile due to the relativistic spin precession of 
pulsar B dBreton et al.ll2008l) . We fitted the high time resolu- 
tion lightcurves to a double Fermi function: 



1 



1 



2 (0+0i)/wj + \ e (.<H<t>,)/w<. + 1 



+ (!-/?), (1) 



where <fi is the orbital phase defined from the superior con- 
junction of pulsar A, and the subscripts i and e refer to the 
ingress and egress, respectively, of the location (</>; ie ) and 
steepness (w,-„) parameters. Since the light curves were nor- 
malized such that the out-of-eclipse values are unity, we 
forced the out-of-eclipse level to be unity and kept the eclipse 
depth, /?, as a free parameter. Hence, /3 = and 1 implies 
full transmission and absorption, respectively. In this frame- 

14 Data at higher frequencies could not be used because they do not have 
sufficient S/N to detect the eclipses. 
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Figure 5. Illustration of the truncated magnetosphere of pulsar B within which the absorbing plasma is confined as seen from an observer located on Earth. The 
horizontal line represents the projected motion of pulsar A, wh i ch goe s from left to right, with the labels indicating the orbital phases in degrees. The system 
configuration is that of the best-fit geometry from Breton et al. 12008), with the spin axis, depicted as a gray rod, tilted 130° with respect to the orbital plane 
(vertical, in the plane of the page), and oriented 51° away from our line-of-sight vector (pointing directly off the page). The angle between the spin and magnetic 
axes is 71°. The four panels display different spin phases of pulsar B, with the upper left corresponding to pulsar B's spin phase 0.0 (i.e., when pulsar B's radio 
emission is illuminating us). Depending on the location of pulsar A, the rotation of the magnetosphere of pulsar B can modulate the flux from A twice per rotation 
at the ingress (leftmost side), or only once closer to the conjunction time (center). 



work, we can evaluate the full width at half-maximum as 
FWHM= <j> e - tpj. Results are presented in Table [3] 

Figure [6] shows the eclipse duration as a function of radio 
frequency inferred from the FWHM for the data reported in 
Table [3] Since the goodness-of-fit is difficult to assess given 
that we fitted a smooth function to profiles that have more 
complex structure, we verified that our inferred errors for the 
duration of the eclipses at 820 MHz are consistent with the 
observed scatter of measured values at this frequencyQ We 
observe that the frequency dependence of the eclipse duration 
is clearly not steep enough to account for FWHM oc i/ -1 / 3 pre- 
dicted by the RG05 model as introduced in Section Q] This 
model relies on the fact that the enhanced absorbing plasma 
density is provided by pulsar A's radio beam illumination and 
cyclotron absorption increases the trapping efficiency of elec- 
trons in the magnetosphere. 

We compared the relationship bet ween the eclipse dura- 
tion an d the radio frequency with the iLvutikov & Thompson! 
(120051) models which are based on synchrotron absorption 
from the relativistic electrons confined within the closed field 
lines of the magnetosphere. In the first version of the model, 
LT05I, a constant electron density distribution is assumed for 
the plasma inside the closed field lines, before dropping to 
zero f or the field lines clo sing beyond the truncation radius 
7? mag . iBreton et alJ (120081) demonstrated the success of this 
model at reproducing the 820 MHz eclipses. Because of the 

15 We chose 820 MHz because it contains the largest number of data points. 



strong dependence of the eclipse profile on the geometry, we 
worke d using the best-fit parameters found in IBreton et alJ 
(2008) and hence the radio frequency was the only parame- 
ter to vary in th e calculation of the ecl i pse du ration. From 
Equation (52) in lLyutikov & ThoTn pson (2005]), we can write 
the optical depth as 



T v = jJLV 



-5/3 



I(x)(B sin k) 2 / 3 dx , 



with 



I(x) = 



1 if Rmaxix) < Rmag and VB, e /u > (je) - 10 

otherwise 



(2) 



(3) 



Here K is the angle between the line-of-sight and the local 
magnetic field direction, /j, is a constant that accounts for var- 
ious physical quantities and plasma properties, R mlix stands 
for the maximum extent of a field line, R mas is the truncation 
radius of the magnetosphere, i/^.e is the local synchrotron fre- 
quency, and (7 e ) = 3ksT e / m e c 2 . 

As we can see in Figure [6j the eclipse duration inferred 
from the LT05I model is able to reproduce the data fairly well, 
though the 1950 MHz observations suggest a longer duration 
than expected. Given the poor S/N and the wide bandwidth 
in this band, it is not clear whether this difference is signif- 
icant or not. In principle, data at higher radio frequencies 
(> 3000 MHz) would provide a good test but the fact that pul- 
sar A becomes very faint beyond 2200 MHz poses an observ- 
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Table 3 

Double Fermi Function Fit Results" 



Date 




4>e 


Wi 


We 


P 


FWHM 


(MJD) 


(deg) 


(deg) 


(deg) 


(deg) 


(Normalized) 


(deg) 


325 MHz 


53191 


0.52(4) 


0.90(0) 


0.22(3) 


0.01(0) 


0.96(3) 


1.42(4) 




0.61(5) 


0.91(1) 


0.29(4) 


0.05(1) 


1.00(4) 


1.53(5) 


427 MHz 


53005 


0.65(6) 


0.83(2) 


0.24(5) 


0.08(2) 


1.00(6) 


1.47(7) 




0.59(3) 


0.85(1) 


0.18(3) 


0.04(1) 


0.97(3) 


1.45(4) 




0.59(4) 


0.85(1) 


0.24(3) 


0.05(1) 


1.00(4) 


1.43(4) 


820 MHz 


52992 


0.55(5) 


0.75(2) 


0.23(3) 


0.07(1) 


0.99(5) 


1.30(5) 




0.58(3) 


0.75(1) 


0.21(3) 


0.05(1) 


1.00(3) 


1.34(4) 


52997 


0.62(3) 


0.76(1) 


0.20(3) 


0.06(1) 


0.96(3) 


1.38(4) 




0.64(3) 


0.76(1) 


0.22(3) 


0.06(1) 


0.99(3) 


1.40(4) 




0.62(3) 


0.73(0) 


0.14(2) 


0.01(0) 


0.91(3) 


1.36(3) 


53211 


0.70(3) 


0.79(2) 


0.10(3) 


0.04(1) 


0.92(3) 


1.49(4) 




0.57(6) 


0.78(2) 


0.19(5) 


0.06(2) 


0.98(6) 


1.35(6) 


53311 


0.61(5) 


0.79(1) 


0.14(4) 


0.01(1) 


0.88(5) 


1.39(5) 




0.56(6) 


0.76(2) 


0.21(5) 


0.06(2) 


1.06(7) 


1.32(6) 




0.54(16) 


0.73(3) 


0.37(10) 


0.07(3) 


1.13(16) 


1.27(16) 


53379 


0.69(5) 


0.79(3) 


0.16(4) 


0.10(3) 


0.96(5) 


1.48(6) 




0.60(6) 


0.76(2) 


0.25(5) 


0.05(2) 


1.02(6) 


1.36(7) 


1400 MHz 


52984 


0.51(7) 


0.76(4) 


0.20(5) 


0.09(3) 


0.91(8) 


1.28(8) 




0.51(5) 


0.69(2) 


0.16(4) 


0.04(1) 


0.94(5) 


1.20(5) 


1950 MHz 


53212 


0.52(7) 


0.67(2) 


0.23(5) 


0.05(2) 


1.01(8) 


1.19(7) 




0.56(3) 


0.64(0) 


0.11(3) 


0.00(0) 


0.87(4) 


1.21(3) 


53378 


0.64(3) 


0.67(0) 


0.13(3) 


0.01(0) 


0.89(3) 


1.31(3) 




0.59(5) 


0.70(2) 


0.18(4) 


0.06(2) 


0.91(6) 


1.30(6) 



"Numbers in parentheses are the lcr uncertainties on the last significant dig- 
its. 

b Since the pulsed flux level is normalized to be unity in the out-of-eclipse 
region, depths /3 of 0.0 and 1 .0 indicate full transmission and total absorption, 
respectively. 

Convergence problems with the non-linear least squares algorithm were 
encountered with the second eclipse of MJD 53379 and hence it was dis- 
carded from the eclipse duration and depth analysis. 



ing challenge. The co-addition of several individual eclipses 
would improve the S/N but the process requires shifting the 
light curves so that t he spi n phases of pulsar B are aligned, as 
described in Section [3.1.11 in order to sum the modulation fea- 
tures coherently. We do not possess sufficient high-frequency 
data to attempt thi s and since pulsar B is now undetectable 
dPerera et al.ll2010l) . it would be technically difficult, if not 
impossible, for eclipses that are not recorded in a continuous 
observation to be combined. A possible workaround, though 
not ideal, would be to conduct another multi-frequency cam- 
paign in the same bands as our current data but including 
much more 1950 MHz observations in the hope of obtaining 
better statistics and perhaps recording some eclipses with low 
radio frequency interference. 

We also examined the possibility of a different distribution 
for both the electron density and temperature. If one assumes 
that the electron density and temperature are constant along 
each field line, and are power-law functions of the distance 



(i.e. 



T e oc i ?°„ T and n R o c K 



A then u sing Equations (39) and 



(50) of ILyutikov & Thompson! (120051) one can write the opti- 
cal depth as 



T V = \XV 



-5/3 



/+oo 
^ m axW^- 5Q/3 )(Bsin K ) 2 / 3 ^, (4) 
oc 
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Figure 6. Eclipse duration as a function of radio frequency for data collected 
between 2003 December and 2005 January. Here, the duration corresponds 
to the full width at half-maximum of a double Fermi function fitted to the 
eclipse profiles (Equation {T}). The vertical error bars represent the lcr un- 
certainties on the fits while the horizontal error bars indicate the receiver's 
bandwidth. The observed behavior of the duration is in good agreement with 
that ca lculated from light curves simulated using thejLyutikov & T hompsor] 
120051) s ynchrotron absorptio n/dipolar geometry model at the best-fit parame- 
ters from Breton et al. 1 2008 ) — LT05I (dashed black) use a constant electron 
density with a hard cutoff radius, while LT05II (dashed red) use an electron 
density that increases like a power law with radius. The spread in the theoret- 
ical predictions of LT05I is due to fact that the modulation features slightly 
change depending on the zero point of the rotational phase of pulsar B with 
respect to the orbital position. The dotted line shows the best-fit i/~'/ 3 rela- 
tionship predicted by the RG05 model. The dash-dotted line shows a hybrid 
version of the LT05I and LT05II models, with a truncated power-law distri- 
bution (see Sectionf4). 

with y! being a constant that accounts for various physical 
quantities and plasma properties. 

One readily recovers Equation (f2]i by assuming constant 
temperature and density distribution (i.e., fixing a = (3 = 0, 
and imposing the truncation radius). On the other hand, the 
LT05II model ( see Section [D proposes that th e electron den- 
sity is given by (ILyutikov & Thompsonl200"5l Equation (56)) 



n e = 



k B T e B 2 
m p c 2 2itm p c 2 ' 



(5) 



w hich implies that (3 = a — 6, w ith a suggested to be equal to 
7 (ILyutikov & Thompsonll2005[ Equation (59)). 

As we can see from Figure [6] the behavior of the LT05II 
eclipse duration is for most parts very similar to that of the 
LT05I and is relatively consistent with the observations. This 
behavior was expected since the only explicit frequency de- 
pendence of LT05I and LT05II comes from the j/ -5 / 3 compo- 
nent of the synchrotron absorption coefficient. 

We note, however, the departure of the LT05II model at low 
frequency, where it produces eclipse durations that increase 
exponentially and are marginally longer than observed. This 
is probably a consequence of the very peaked electron power- 
law temperature distribution, which produces a runaway ab- 
sorption of low-frequency radio photons at large magneto- 
spheric radii. Finally, we also find that the eclipse duration 
is not very sensitive to a variation of the power-law index. 
This is mainly due to the constant -6 factor arising from the 
B 2 dependence in the electron density (see Equation (0) of 
the LT05II model, which diminishes the importance of the 
power-law index a. 

3.3. Depth 

Complementary to the eclipse duration in the investigation 
of the plasma properties is the radio frequency dependence of 
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Figure 7. Same as in Figure [6] but for the eclipse depth as a function of 
radio frequency for data collected between 2003 December and 2005 January. 
The depth corresponds to the parameter /3 in Equation (TV which fits the 
eclipse light curve for a double Fermi function (see TableQ). The vertical 
error bars are the \a uncertainties on the fits and the horizontal error bars 
indicate the receiver's bandwidth. The dashed lines represent the prediction 
of two versions of the Lvutikov & Thompson 1 2005) model, with a constant 
plasma density and sharp cutoff radius in black, LT05I, and with a radial 
power-law distribution of the electron density and temperature in red, LT05II. 
The dash-dotted line shows a hybrid version of the LT05I and LT05II models, 
with a truncated power-law distribution (see Sectionf4). 

the eclipse depth, which we define as the j3 parameter in Equa- 
tion (QJ. Figure [7] suggests that there is little or no change 
in the eclipse depth as a function of frequency in the range 
325- 1950 MHz. This contrasts with the expected behavior 
of the hard cutoff model of LT05I which predicts a rough lin- 
ear decrease of the depth between 1000 and 3500 MHz. Con- 
versely, the power-law distribution of the LT05II model dis- 
plays much better agreement as its radio frequency evolution 
is very mild. 

4. A HYBRID MODEL 

The LT05I and LT05II models both display contradictory 
features. On the one hand, LT05I agree fairly well with the 
duration of the eclipses, whereas LT05II predict durations that 
are too long in the low radio frequency regime. On the other 
hand, LT05II seem to account much better than LT05I for the 
very small variation of the eclipse depth as a function of fre- 
quency. The main conclusion that should be drawn is that the 
magnetosphere must possess a region where the optical depth 
is large over a wide range of radio frequencies before drop- 
ping abruptly in a very sharp transition region. It therefore 
appears that the actual distribution of temperature and density 
might be a compromise between the LT05I and LT05II mod- 
els: a power-law which is truncated at a characteristic radius, 
i.e. T e oc n e oc r a for R max (r) < R mag and otherwise. In this 
fashion, the optical depth would be sufficiently large to high- 
frequency radio photons traveling within the cutoff radius to 
produce deep eclipses even in this regime. At the same time, 
the truncation of the magnetosphere would limit the spatial 
extent of the eclipses at low frequencies where electrons have 
larger synchrotron absorption cross-sections to low-frequency 
photons over greater radii. Figures[6]and[7]show such a hybrid 
model. Without any attempt to optimize the power-law index 
and the truncation radius, it displays much better agreement 
to the data than the other models. 

We note that this hybrid model shares similarities with 
Earth's magnetosphere which is immersed in the solar wind 
and possesses a sharp transition region, the plasmapause, 
where the plasma densi ty drops by almost two orders of mag- 
nitude (Carpent erl 19631) . The sudden decrease of plasma den- 



sity is due to the fact that confinement is efficient only where 
the field lines are near corotation. Instead of extending all the 
way to the limit of the light cylinder, as one would naively 
expect for an isolated radio pulsar, the plasmasphere of pul- 
sar B is confined closer to the neutr on star. As demons trated 
by the particle-in-cell simulations of lArons et al.l (120051) . pul- 
sar A's magnetized wind penetrates deep into the magneto- 
sphere of pulsar B and reconnects with its magnetic field. 
This generates a tangential stress that breaks down corota- 
tion. We can compare the size of the plasmasphere that we 
obtain from our modeling with that estimated from the stand- 
off distance determined by equating pulsar A's wind pres- 
sure to pulsar B's magnetic press ure. The inferred plasma- 
spher e radius from our modeling dBreton et a l. 2008; Bretonl 
120091) is R mRg = 0.023 Dab, where Dab is the orbital separa- 
tion. This implies a physical size of 2 x 10 9 cm, which is 
half that estimated using the standard magnetic di pole braking 
equat i on with fiducial neutron star parameters (lArons et alj 
120051: iLvutikov & Thompsonll200"5l) . Such a conundrum is 
not impossible to resolve if pulsar B/A's magnetic field has 
a smaller/larger strength than estimated from the spin-down, 
and/or if further numerical simulations can show that the spin- 
down torque on pulsar B can be larger than that inferred for a 
dipole in vacuum. If arising from the torque only, this would 
require a factor ~ 16, which is not incompatible wit h the limit 
of ~ 10 estimated by lLyutikov & Thompson! (120051) . 

From our modeling of the eclipses we observe that a 
rather large plasma density is required in pulsar B's magne- 
tosphere, with an electron multiplicity factor ~ 10 5 the stan- 
dard Goldreich-Julian density. One possibility to explain the 
high plasma density is that it is specific to the double pul- 
sar system as a consequence of the interaction between pul- 
sar B's magnetosphere and pulsar A's relativistic wind. "Ra- 
dial diffusion" type effects such as those observed in plane- 
tary magnetospheres could bring a high particle influx insid e 
the magnetosphere of pulsar B (ISchulz & Lanzerottilll974l) . 
On the other hand, another possibility might be that a large 
electron multiplicity is instead a generic feature to all pulsars 
(e.g., using diffusion from the open field lines, for example; 
lArons & Smithl[T979l) . More detailed theoretical would be re- 
quired to study the viability of different mechanisms at repro- 
ducing large plasma density. However, without secure mea- 
surements for other sources it is difficult to draw conclusions 
for the pulsar population as a whole. 

5. CONCLUSION 

In this paper, we have presented the detection of flux mod- 
ulations in the eclipse light curves of pulsar A at radio fre- 
quencies betwee n 325 and 1950 MHz, exte nding from the ini- 
tial discovery of lMcLaughlin et al.l d2004al) at 820 MHz. This 
shows that this behavior is intrinsic to the eclipse mechanism. 
Furthermore, our quantitative analysis of the light curves re- 
veals that the periodicity is indeed exactly related to the period 
and higher harmonics of pulsar B's rotation and only occurs 
within the orbital range of the eclipse. The flux modulations 
show a dynamic evolution, switching between modes domi- 
nated by power at the spin period of pulsar B and that of its 
first harmonic. Their phases and behavior are consistent with 
the absorbing material corotating with pulsar B and a dipolar 
magnetic field having a sharp optical depth transition region 
located in the closed field lines. 

Regarding the prope rties of the absorbing plasm a, we find 
that the RG05 model dRafikov & GoldreichlEOOBT) produces 
an eclipse duration that is too radio-frequency sensitive. This 
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disfavors the hypothesis that a only "modest" electron multi- 
plicity (~ 100) is reached from the illumination of pulsar A's 
radio flux and sustained in a steady state by the heating via 
resonant cyclotron absorption. Instead, our observations sug- 
gest that the plasma contained in pulsar B's magnetosphere 
possesses a "large" electron multiplicity (~ 10 5 ). 

Unfortunately, the density and temperature distributions of 
the relativistic electrons are more difficult to constrain. From 
our work, we conclude that for a wide range of radio pho- 
tons the optical depth is very large up to a characteristic cut- 
off radius beyond which it drops abruptly. Such a behav- 
ior would be reproduced by a hybrid model of LT05I and 
LT05II for which the plasma electron temperature and den- 
sity distributions are both truncated radial pow er-laws and 
resemble Earth's plasmasphere as suggested by lArons et all 
( 120051) . The measured size of the plasmasphere 2 x 10 9 cm is 
about half that of the expected magnetic field to wind pressure 
standoff distance and could imply that the spin-down of pul- 
sar B is larger than expected from standard estimates. It is not 
clear whether the properties of pulsar B's magnetosphere can 
be generalized to other isolated pulsars since interactions with 
pulsar A's wind are certainly not negligible. Nevertheless, it 
represents one of the best available opportunity to study the 
immediate vicinity of a pulsar. 
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